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■ Abstract 

O 

o. 

CN . We construct a two parameter family of collapsing solutions to the 4+1 Yang- 

Mills equations and derive the dynamical law of the collapse. Our arguments 
indicate that this family of solutions is stable. The latter fact is also supported 
by numerical simulations. 

1. Introduction 



> 

Blow up problems for nonlinear Schrodinger, wave and heat equations have been a 
subject of active research in the last 15 years (see [1,2,3] for reviews and [4,5] for recent 
papers on the subject) . A further surge of interest in blow-up for nonlinear wave equations 
has been recently motivated by their role in attempting to understand the problem of 
singularity formation in General Relativity (see [6] for a recent review). In this paper we 
describe the asymptotic dynamics of blowup for radial solutions of the semilinear wave 
equation 

> ■ 

S3 1 i 

u = Au+-f(u) 1 (1.1) 

^ ■ in iR 2 , where u = u(t, r), r is the radial variable and 

/(«) = 2u(l - u 2 ). (1.2) 

Our analysis is applicable to a wider class of "double- well" type of nonlinearities f(u) pro- 
ducing kink-type solutions, though some of such nonlinearities, e.g. wave maps nonlinearity 
f(u) = — \ sin(2it), lead to certain subtleties and will be considered elsewhere. 

Before stating the results we show how the equation (1.1) arises and put the prob- 
lem in a broader context. We consider Yang-Mills (YM) fields in d + 1 dimensional 
Minkowski spacetime (in the following Latin and Greek indices take the values 1, 2, . . . , d 
and 0, 1, 2, . . . , d respectively). The gauge potential A a is a one-form with values in the 
Lie algebra g of a compact Lie group G. In terms of the curvature F a p = d a Ap — dpA a + 
[A a , Ap] the YM equations take the form 

d a F afi + [A a ,F^] = 0, (1.3) 
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where [, ] is the Lie bracket on G. For simplicity, we take here G = SO(d) so the elements 
of g = so(d) can be considered as skew-symmetric d x d matrices and the Lie bracket is 
the usual commutator. Assuming the spherically symmetric ansatz [7] 

Aii(x) = (5^ - 6^) 1 ~ U ^ r \ (1-4) 

equations (1.3) reduce to the scalar semilinear wave equation for the magnetic potential 

u(t, r) 

u = A (d _ 2) w + ^—^-u{\ - u 2 ), (1.5) 

where A( d _ 2 ) = d 2 + ^-^-d r is the radial Laplacian in d — 2 dimensions. 

The central question for equation (1.5) is: can solutions starting from smooth initial 

data 

"(0,r) = /(r), u(0, r) = g(r) (1.6) 

become singular in future? An answer to this question depends critically on the dimension 
d. To see why we recall two basic facts. The first fact is the conservation of (positive 
definite) energy 

oo 

E = j ' (ii 2 + u' 2 + ^(1 - u 2 ) 2 ^ r d ~ 3 dr. (1.7) 
o 

The second fact is scale invariance of the YM equations: if A a (x) is a solution of (1.3), 
so is A a (x) = \~ 1 A ce (x/\), or equivalently, if u(t,r) is a solution of (1.5), so is u(t, r) = 
u(t/X,r/X). Under this scaling the energy scales as E = \ d ~ A E : hence the YM equations 
are subcritical for d < 3, critical for d = 4, and supercritical for d > 5. In the subcritical 
case, shrinking of solutions to arbitrarily small scales costs infinite amount of energy, so it is 
forbidden by energy conservation. This is a heuristic explanation of global regularity of the 
YM equations in the physical dimension which was proved in [8] and [9]. In contrast, in the 
supercritical case shrinking of solutions might be energetically favorable and consequently 
singularities are anticipated. In fact, for d > 5 equation (1.5) admits self-similar solutions 
which are explicit examples of singularities [10,11] and numerical simulations indicate that 
the stable self-similar solution determines the universal asymptotics of blowup for large 
initial data [12]. 

In the critical dimension d = 4 the problem of singularity formation is more subtle 
because the scaling argument is inconclusive. In this case there are no smooth self-similar 
solutions, however there is a family of static solutions x(r/A), where A > and 

X(r) = (1-8) 
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Using physicists' terminology we shall refer to this solution as the instanton. Numerical 
simulations indicate that the existence of the scale-free instanton plays a key role in the 
dynamics of blowup, namely the blowup has the universal profile of the instanton which 
shrinks adiabatically to zero size [12]. More precisely, it was conjectured in [12] that near 
the blowup time t* the solution has the form 



where the scaling parameter X(t) tends to zero as t — > £*. A natural question is: what 
determines the evolution of the scaling parameter; in particular, what is the asymptotic 
behaviour of X(t) for t — > £*? In this paper we address this question and show that for 
some initial conditions that are close to the instanton 



as t — > t*. The logarithmic correction to the self-similar behaviour is characteristic for the 
blow-up in critical equations - it implies that the speed of blow-up goes asymptotically 
to zero and consequently no kinetic energy concentrates at the singularity (for a different 
approach see [13]). 

2. Results 

Thus we consider the intial value problem (1.1), (1.2) and (1.6). Since we consider 
radial solutions only the full Laplacian A can be replaced by the radial Laplacian A r = 



Our main result states that if initial conditions (1.6) are sufficiently close to 




(1.9) 




(1.10) 



\d T rd r . 




where Aq > and Ao < 0, then the resulting solution is of 



the form 




(2.1) 



where the scaling parameter A = X(t) satisfies the following equation 



AA=?A<, 



(2.2) 



with the initial conditions Ao and Ao . In fact, our procedure allows us to find the solution 
u(x, t) to any order in A 2 with the term of order A 2 given explicitly. 



Note that solutions of Eqn (2.2) with the initial conditions such that Ao < 0, decrease 
to zero as t — > t* for some £* with |A| decreasing so that our approximation improves as 
t — > t*. This and Eqn (2.1) imply that the instanton collapses as t — > £*. 
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To demonstrate the property of Eqn (2.2) mentioned above we note that Eqn (2.2) 
can be integrated explicitly. Indeed, rewrite (2.2) as A _3 A = §A _1 A and integrate the 
resulting equation to obtain 

-lA- 2 = ^(lnA + lnc), (2.3) 



where c > 0. The latter equation can be rewritten as 



A2 MsK (2 - 4) 



This relation shows that we must have 

cA < 1. 

Using Eqn (2.4) we obtain the equation for c in terms of A 



We have two cases. 



a. A > 0. Then cA f 1 and A | oo as t approaches some t* > 0. Moreover, c A = 1 — 
(§c 2 ) 1/3 (**-*) 2/3 as*-**. 

b. Ao < 0. Then A < for t > and there is £* > s.t. A — as t | £*. The value t* 
can be found from (2.4): 

= /l/ rfAlnl/2 (a) ■ < 2 ' 6 > 

o 

Taking into account (2.5) this gives 

^wAolAol" 1 . (2.7) 

The time is the point of collapse. 

1/2 

Note that in this case the function |A| decreases as t — £* as [in (^)j and there- 
fore our approximation improves as t — t*. 

Solutions of Eqn (2.2) with Ao < and Ao > have the following asymptotics as 

X = ./| «•-« . (2.8) 
V 3 ^- ln(i, - i) 

In conclusion we observe that Eqn (2.2) is invariant under the transformation 

X(t) — /j~ 1 X(/jt) 

inherited from the invariance of parent Eqn (1.1) under the scaling transformation 

u(r, t) — u(fj,r, fit). 
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3. Scaling transform and zero mode 

A key role in our derivation is played by the fact that Eqn (1.1) is scale covariant 
under the transformation 

u(r,t) ->u(r/A,t/A), (3.1) 

i.e. if u(r, t) is a solution to (1.1), then so is u(r/X, t/X). In particular, if v(r) is a stationary 
solution, then so is v(r/\), A > 0. The infinitesimal change of the instanton x under this 
transformation is 

X - X + (3.2) 

where the function ( is defined by 

C(r) := d x \x=ix(r/X) = -rd r x(r). (3.3) 

Explicitly 

CM = (1^)1 ■ (3-4) 

Of course, ( is the zero mode, 

L( = 0, (3.5) 
of the linearization of the r.h.s. of (1.1) on x, i.e. of the operator (recall, A r = ^d r rd r ) 

L:=-A r -l/'( X (r)). (3.6) 

(This operator is the variational or Frechet derivative, L = <90(x), of the map <p{u) = 
—Au — -pzf(u) at the instanton %.) 

The following properties of the operator L will be important for us: 

- L = L* > 

- L has a simple eigenvalue at with the eigenfunction £ 

- the continuous spectrum of L fills [0, oo). 

The first and third properties are obvious and the second property follows from the 
equation L( = and the fact that ( > by the Perron- Frobenious Theory (see [14,15]). 

Consider solutions of Eqn (1.1) of the form u(r, t) := v(r/X,t), where A > depends 
on t. Plugging the function u(r, t) = v(r/X,t) into Eqn (1.1), we obtain the following 
equation for v and A: 

A y v + y~ 2 f(v) = -X 2 B lV - XXB 2 v + X 2 dfv - 2XXB 2 d t v, (3.7) 
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where 

Bi = -yd y - (yd y ) 2 and B 2 = yd y . (3.8) 

4. Orthogonal decomposition 

We look for a solution of Eqn (1.1) of the form 

u(r, t) = v(r/X, t) = x(r/A) + w(r/X, t) (4.1) 

with A = X(t) and 

w(y,t) small for all times. (4-2) 

Moreover, to fix the splitting between the dynamics of A and of w we require that w is 
orthogonal to the zero mode (: 

oo 

J ay)w(y,t)ydy=0. (4.3) 
o 

The last two conditions will give us the dynamic law for A. 

Now we plug the decomposition v(y,t) = x(u) + w (y^) into Eqn (3.7) and use that 
the function x satisfies the equation 

A yX + y- 2 f(x) = 0, (4.4) 

to obtain the equation for w: 

(L + X 2 d 2 )w = F(w,X), (4.5) 
where, recall, L = L x is the linearized operator around x 

L:=-A y -y- 2 f(x) (4-6) 

and 

F(w, A) := X 2 B 1 ( X + w) + XXB 2 ( X + w) 

+ y- 2 N(w) + 2XXB 2 d t w (4.7) 

with the nonlinearity N(w) defined by 

N(w) := f(x + w)- f(x) ~ f(x)w, (4-8) 
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which in the case f(u) = 2(1 — u 2 )u gives 

N(w) = -6xw 2 - 2w 3 . (4.9) 

5. Perturbative analysis. Outline 

We explain the main idea of our approach by proceeding formally and ignoring infrared 
divergences arising in an attempt to justify our analysis. In the next section we present a 
full perturbation theory. We look for a solution to Eqn (4.5) in the form 

w(y^) = J2^ J ^(y)- (5-1) 

i>i 

Plugging this expansion into Eqn(4.5) we arrive at a series of equations 

L^ = F i (&,...,^-i), (5-2) 

j > 1, where £ = X- 

We demonstrate our approach by analyzing the cases j = 1 and 2 in detail. We begin 
with j = 1. It is clear from (4.5)-(4.7) that 

F 1 (£o) = B lX . (5.3) 

Thus £i satisfies the equation 

L£i = B lX . (5.4) 

Since, as we show in Appendix 1, 

f (Bi X = 0, (5.5) 
Eqn (5.4) has a solution. The general solution of this equation is 

£i(y) = £io0/) + aiC(2/) (5-6) 

for any ol\ G M. Here 

Cio(y) = - (1 ^ 2)2 (5-7) 

Now plugging m; = A 2 ^i + 0(A 4 ) into (4.3), and using (5.5) and B 2 x = we obtain 
J C (A 4 Si6 - AAC - A 4 6y- 2 X ei 2 ) = 0(A 6 ) + 0(AAA 2 ). (5.8) 
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It will be shown in Appendix 1 that the coefficients in front of a\ and a\ (remember (5.6)) 
vanish: 

J CBiC - 12 f Cy" 2 x6oC = (5.9) 

and 

/ Cy-\C 2 = o (5.io) 

Therefore relation (5.8) becomes 

AA- 7 A 4 = 0(A 6 ), (5.11) 

where 

7 = J C (Bi£io " 6y- 2 X^ 2 o) (| C 2 ) 1 • (5-12) 
We show in Appendix 1 that 7 = 3/4 which, in the leading order, brings us to Eqn(2.8): 

AA = ^A 4 . (5.13) 

As (5.11) shows, this equation is valid modulo the correction 0(A 6 ). 

Now we proceed to the second term in (5.1) and derive a correction to Eqn(5.11) 
(or(5.13)). Remember that £2 is defined by (5.2) with j = 2. Keeping in mind Eqn(5.13) 
we have 

F 2 (£ ,£i) = B 1 Z 1 + 1 B 2X -6y- 2 X ei- (5-14) 

By (5.8) we have that 

J C^CcCi) = (5.15) 

so that the equation L£ 2 = F 2 (^o^i) is solvable. Now plugging 

w = A 2 ^ + A% + 0(A 6 ) (5.16) 

into (4.3), we obtain the equation for A 

AA- 7 A 4 + 5A 6 +eA 2 AA + 0(A 8 ) =0 (5.17) 

with 6 and e given in terms of integrals of £1 and £2 which can be explicitely computed. 
Observe that Eqn(5.17) is equivalent to the equation 

AA - 7A 4 + (5 - 7 e)A 6 + 0(A 8 ) =0. (5.18) 
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We can continue in the same manner to find the equation for A to an arbitrary order in 
A 2 . 

Though the perturbation theory outlined above leads (as we will see in the next 
section) to correct — in the leading order — equations for the dilation parameter A, it is, in 
fact, inconsistent. The leading correction, £1, does not vanish at infinity and consequently 
the resulting solution has infinite energy. Worse, higher-order corrections grow at infinity. 
Moreover, orthogonality condition (4.3) is not applicable (and as a result the parameter 
cki in (5.6) cannot be determined). The reason for this inconsistency is that the term 
—X 2 dfW cannot be treated as a perturbation at large distances. A correct perturbation 
theory taking into account the leading contribution of this term at infinity is presented in 
the next section. 

6. Perturbative analysis 

In this section we justify formal analysis of Section 5. We look for a solution, w, of 
Eqn (4.5) in the form 

w(y,t) = J2^ 2j Cj(y)^ 4 y 2 ,t). (6.i) 

i>i 

We fix the functions £j and cpj by requiring that (a) £j and cpj are of the order 0(1), (b) 
the functions cpj satisfy the relations 

<Pj(z,t) = l for z<l (6.2) 

and 

<p j (z,t) = <p j (z)+0(\ 2 ), (6.3) 

(c) the following equations are satisfied 

(L + X 2 d 2 )(X 2 ^ m ) = X^F j , (6.4) 

where the functions Fj are 0(1) and depend only on ^o^o, • • • ,£j-i<fij-i with £o = X an d 
ip = 1: 

F j = FjitoVo, ■ ■ .,€j-i<pj-i) (6.5) 
and (d) the functions £j(y) satisfy the equations 

L^ = F i (&,...,^-i). (6-6) 

As will be shown below these requirements will define the functions £j and <pj uniquely, at 
least in the leading order. 
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We demonstrate our approach by analyzing the cases j — 1 and 2 in detail. We begin 
with j = 1. It is clear from (4.5)-(4.7) that 

Fifow) = B 1X . (6.7) 

Thus £i satisfies the equation (5.4). Recall that due to (5.5) the latter equation is solvable 
and its general solution is given by (5.6). The constant ol\ in (5.6) is determined from the 
condition 

V£m = 0. (6.8) 

Since it plays no role in what follows we do not compute it here (see, however, (6.23) 
below) . 

Now plugging w = A 2 £i</?i + (3(A 4 ) into (4.3), omitting (pi (justification for this will be 
provided later) and using (5.5) and B 2 x = — C we obtain (5.8) which as shown in Section 5 
leads to (5.11) with 7 = 3/4. 

Now we return to expansion (6.1) and find the equation for cpi. Recall that <pi is 
defined through equations (6.2)-(6.6) with j = 1. We derive from these equations the 
equation for (pi(z, t) in the leading order in A 2 and in the domain y 1. To this end we 
use Eqn (5.13) to estimate the order of higher derivatives of A. In the leading order we 
can ignore the dependence of <pi(z, t) on t. Using Eqns (6.4) with j = 1 and Eqn (6.7) and 
using that for y 3> 1 

£i = -l + o(^), yd y ^ = -^ + 0^y B lX = -^+0^y (6.9) 
we obtain after lengthy but elementary computations that cpi satisfies the equation 



Z 2 d\ip X + (z + 7Z 2 ) d z ip! - (l - (f! = -1. 



(6.10) 



To this we add the boundary conditions 

V?i(0) = l. (6.11) 

The second boundary condition, <fii(0), or, alternatively, an arbitrary constant in the 
general solution to (6.10)-(6.11), is found by matching the solution to (6.10)-(6.11) in the 
region z C 1 with solution to (6.4) with j = 1 and (6.7) in the region y ^> 1. This is done 
in Appendix 2 where it is also shown that 



(lnA-l) forz«l 
Vl (z) ={ V } (6.12) 



1 _ 2£ 

4 



* + # + *e- 72 for,»l 
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for some constants c and c. Thus we have 

A 2 6(2/Vi(AV) = 0(y- 1 ) for X 2 y » 1. (6.13) 
This implies, in particular, that the integral in (6.8) converges and it gives 

ai = o(lni). (6.14) 

Now we proceed to the second term in (6.1) and derive a correction to Eqn(5.11) 

(or(5.13)). Remember that £2 is defined by (6.6) with j = 2. Keeping in mind Eqn(5.13) 
we choose 

F2{^o,^i) = B 1 {^ 1 ) + ^B 2X -6y- 2 X {^i) 2 . (6.15) 
By (5.8) we have that 

' CF 2 (Zo,ti) = (6.16) 



so that the equation = ^2(£o>£i) is solvable. Eqns(6.2)-(6.6) with j = 2 imply an 
equation for cp2 which is analyzed in a similar way as the equation for cpi. Now plugging 

w = \ 2 £ 1 <p 1 + \ 4 &<p 2 + 0(\ 6 ) (6.17) 

into (4.3) and setting ip\ and <^ 2 to 1, we obtain the equation (5.17) (or (5.18)) for A. We 
can continue in the same manner to find the equation for A to an arbitrary order in A 2 . 

Conclusion 

We found, perturbatively, a two parameter family of collapsing solutions (parametrized 
as Ao and Ao) to the nonlinear wave equation (1.1)— (1.2) arising from the Yang-Mills 
equation in 4+1 dimensions. We also found the corresponding dynamics of collapse. The 
perturbation theory developed suggests that this family is (asympototically) stable. This 
conclusion is supported by numerical simulations we performed (some of the results of 
these simulations are given in the figure below). 
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Appendix 1 Computation of various integrals 

In this appendix we show (5.9), (5.14a), (5.14b) and 7 = 3/4 (see Eqn (5.17)). 
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1. (5.9). Recall C = -yd y x(y) and B x = -yd y - (yd y ) 2 . Hence B x \ = (1 + yd y )( and 



J (B lX = J aC + vC) = J vCivO'dy. 

o 

Integrating by parts we get 



oo 

_. oo 

y((y(Ydy= -{yO 2 











and since the boundary term vanishes, we get J C-^iX = 0- Nota bene, this shows that 
the orthogonality condition for j = 1 is basically equivalent to the square integrability 
of C- 

In what follows we use the following relation (m < n — 2) 

oo 

x m dx m{m — 1) • ... • 1 



(l + x) n (n- l)(n-2) • ... • (n-m- 1) 



2. (5.14a). Show that 

J (B 1 ( -12 J Cy- 2 XWl C = 0. 

Compute 

oo 

CBiC = - J(tt'y 2 + (d y (yC)y 2 dy 
o 

oo 

= /(ccV + cV)^ 



oo 

= / (-C 2 y + C'V 

o 

! / o f ( -y* , V(i- y 2 ) 2 ^ 



This gives 



/cftC=2/( 



(1 + !, 2 ) 4 (l + </ 2 ) 6 



oo 

2 



-/( 



x 2 4x 2 (l - xY . 



(1 + x) 4 (l + x) ( 



2-1 4-2-1 8-3-2-1 4-4-3-2-1 1 
+ „ : = - z : „ „ + 



3-2-1 5-4-3 5-4-3-2 5-4-3-2-1 5' 
12 



BOS - February 7, 2008 





oo 



J (l + x) 7 



3-2-1 4-3-2-1 1 
+ 



6-5-4-3 6-5-4-3-2 60 

Hence J (B 1 ( - 12 J Crt^iC = 8(^ - ^) = 0. 
3. (5.14b). Compute 



43/ 



00 

4 1 „,2 





OO 

x 2 (l — x) 
(l + x) 7 



2-1 3-2-1 



6-5-4 6-5-4-3 



0. 



4. 7 = 3/4. Compute 



00 

1 /V2 o f y 4 



C = 2 / 77^4^ 





00 



-/ 



x 2 dx 

(l + x) 4 


2-1 1 



3-2-1 3 



Compute 



so 7 = f . 



J (B lWl - 6 y (y- 2 xwl = 2, 



Appendix 2 Solution <p-\ 



In this appendix we find the solution to 'initial' value problem (5.19)-(5.20) matching 
the solution to (5.4) with j = 1 in the region 1 < y C A -2 . In the region {z <C 1} Eqns 
(5.19)-(5.20) have the general solution 

<P! = 1 - I lnz(z + . . .) + c'^z(l + ...) (A2.1) 

with an arbitrary constant c'. 



13 



BOS - February 7, 2008 

For z»l Eqn (5.19) has the general solution 

2 1 e~ 72 
= — + + (^2.2) 

7-2 v- 2 V- 2 

with arbitrary constants c and c. Here and are solutions of the corresponding 
homogeneous equation in the region z >> 1. 

It remains to find the constant d in (A2.1). To this end we match (pi(X 4 y 2 ) (in the 
leading order) to the solution of the equation 

(l + \ 2 ^)w = \ 2 B iX (A2.3) 

in the region 1 C <C A -2 . We find the solution of the latter in the leading order in A 2 
by a pertubation theory: 

w = X 2 w 1 +X 6 w 2 + (A2A) 
where w\ = £i (see Eqn(5.6)). This implies the equation for w 2 : 

Lw 2 = 2 7 £ 10 . (A2.5) 
Two solutions of the corresponding homogeneous equation are (see (3.5)) 

y 2 _ y 2 3 13 1 3y 2 lny 2 

^"(TTTF m ~ T + 2 " 4(PTT) " V(y 2 + 1) + W+W { } 

By the method of variation of constants we obtain 

w 2 = c\7]\ + c 2 i] 2 (A2.7) 
where the functions c\ and c 2 are given by 

y\ 2 4 1 y 6 lny 2 



a = -is + y ~ 9 , ; + 



y 2 + l (y 2 + l) 2 (y 2 + l) 3 



3y 4 lny 2 3y 2 lny 



2(y 2 + l) 2 y 2 + 



9 ^ 

y f ds In s , A „ 

r + 3 /7TT + 4 (A2 ' 8 > 

n •/ 



and 



C2 = 7 < ln(/ + 1) + ^__^ + _l__iI>. (J 4 2 . 9 ) 
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Eqns (A2.7)-(A2.9) for y > 1 yield 



w 2 = -|?/ 2 + 7^7 ln (! + 2/ 2 ) + • • • 
= |y 2 lny 2 + ... 



(A2.10) 



Since on the other hand wi = —1 + 0(y 2 ), we find in y ^> 1 that 



«; = -A 2 



1 — Inz + 
4 



(A2.ll) 



where z = X 4 y 2 . Comparing (A2.ll) with (A2.1) we find 



c' = Iny 2 -ln(AV) =lnA 



(A2.12) 



and therefore 



pi(z) = 1 - ^zln4- + ... for 
v ; 4 A 4 



(A2.13) 
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0.35 




0.1 - 



0.05 - 



j i i i i i i i i i i 

-60 -55 -50 -45 -40 -35 -30 -25 -20 -15 -10 

In(f-t) 



Comparison of the numerically computed scaling parameter divided by t* — t with the 

analytic formula = ln(£* — t)) -1 / 2 . Note that there are no free parameters 

to be fitted. We believe that the apparent discrepancy (which is of the order of 10% at 
ln(t* —t) = —60) can be accounted for by including higher order corrections to the formula 
(1.10). ' 
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